The Origin of Cosmic Dust. 



m 
o 
o 



Loretta Dunne*, Stephen Eales*, Rob Ivison^, Haley Morgan*, Mike Edmunds* 

* Department of Physics & Astronomy, Cardiff University, 5 The Parade, Cardiff CF 24 3YB, 

UK 

^Astronomy Technology Centre, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 

Large amounts of dust (> lO'^M©) have recently been discovered in high redshifts 
' quasars^'^ and galaxies^^, corresponding to a time when the Universe was less than 
one-tenth of its present age. The stellar winds produced by stars in the late stages of 
^ \ their evolution (on the asymptotic giant branch of the Hertzsprung- Russell diagram) are 
though to be the main source of dust in galaxies, but they cannot produce that dust on 
^ \ a short-enough timescale^ (< 1 Gyr) to explain the results in the high-redshift galaxies. 
^ ■ Supernova explosions of massive stars (type II) are also a potential source, with models 
^ ! predicting 0.2 4 M0 of dust'^^'^. As massive stars evolve rapidly, on timescales of a few 
2 ■ Myr, these supernovae could be responsible for the high-redshift dust. Observations^^~^^ 
r~| . of supernova remnants in the Milky Way, however, have hitherto revealed only 10"'^ — 

o 

' Cassiopeia A. This observation implies that supernovae are at least as important as 



10 ^ Mq of dust each, which is insufficient to explain the high-redshift data. Here we 
report the detection of ~ 2 — 4Mq of cold dust in the youngest known Galactic remnant, 



stellar winds in producing dust in our Galaxy and would have been the dominant 



H ' source of dust at high redshifts. 

Over the past three decades, many searches for dust in supernova remnants (SNR) have been 
made in the mid and far-infrared (6-lOO^m). Remnants must be studied when they are young, 
before they have swept up large masses of interstellar material which makes it difficult to distinguish 
dust formed in the ejecta from that present in the ISM prior to the explosion. The handful of 
Galactic remnants which are both young and close enough (Gas A, Kepler and Tycho) have been 
studied with the Infrared Astronomical Satellite (IRAS) and the Infrared Space Observatory (ISO), 
but although dust at 100-200 K has been detected, the dust mass deduced is only 10""^ - IQ-^Mq, 
many orders of magnitude lower than the solar mass quantities predicted. ^^^^^ The formation of dust 
in the recent supernova 1987A has been implied indirectly from the fading of the silicate line and 
increase in 10/im emission, although the quantity is heavily dependant on the assumptions made 
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about the dumpiness of the dust*. The IRAS/ISO observations were not sensitive to cold dust at 

< 25 K, if the models are correct about the amount of dust produced in supernovae (SNe) then the 
bulk of the dust must be cold. Such dust will emit most of its radiation at longer wavelengths and 
so the best place to search for cold dust emission is at sub-millimetre (submm) wavelengths (0.3-1 
mm) . Here, we present our analysis of submm data for Cassiopeia A, the youngest known Galactic 
remnant. Cas A is the brightest radio source in the sky, and still produces significant synchrotron 
emission at submm wavelengths. Cas A is believed to be the remnant of the explosion of a massive 
(20-30 Mq) progenitor star, which occurred around 320-340 yrs ago at a distance of 3.4 kpc.^^'^^ 

Cas A was observed with the SCUBA^^ array on the JCMT in June 1998. The 850/xm image is 
presented in Fig. 1 and shows a ring-like morphology, similar to the X-ray and radio images-^^'^^. 
Over two-thirds of the 850/im emission in the main ring is synchrotron produced by rclativistic 
electrons spiralling in the intense magnetic field. Fig. 2 presents the spectral energy distribution 
(SED) of Cas A from the radio to the mid-IR, showing the clear excess due to dust at wavelengths of 

< 850)um. The synchrotron radiation displays a constant power-law behaviour over more than two 
decades in frequency with S^, oc where a = —0.72, meaning that its contribution to the submm 
flux can be easily estimated and subtracted. The 450/xm image is shown in Fig. 3, the synchrotron 
contribution is now only one third and emission from cold dust dominates which is why the 450 
and 850^m images appear different. Wc can remove the synchrotron contribution from the submm 
maps by scaling an 83 GHz image, which is the closest in frequency and resolution to our images. 
Fig. 4 shows the result with contours denoting the 450/im synchrotron subtracted flux. The general 
similarity between the two wavelengths, once the synchrotron has been removed, is strong evidence 
that this is indeed emission from cold dust. 

The rings on Fig. 4 indicate the position of the forward and reverse shocks as determined from 
Chandra X-ray data.^^ Most of the dust appears to be contained between the two, where the gas 
density as traced by the X-rays is greatest. We fitted a two-temperature grey body SED to the 
synchrotron corrected IR/submm fiuxes (Fig. 2) determining temperatures of 112 K and 18 K for the 
two components. The SED parameters and their uncertainties are given in Table 1. The temperature 
of the hot dust is easily explained if it is co-extensive with the X-ray emitting gas (10^ — 10^ K) 
and is heated by collisions with fast moving electrons and ions^^. The cold dust could be explained 
in a number of ways and a more detailed investigation will be presented elsewhere (Dunne et al. 
in prep). If the grains are very small (< 0.005/xm) they may cool to 15-20 K in the time between 
collisions with the gas particles. However, the dust SED might then be expected to show dust at all 

*If the dust is uniform, the observations imply a very low efficiency of condensation of heavy elements into dust 
grains (~ 0.001%), while if dust is very clumped this could rise to almost 100% in some species 
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temperatures between 100-20 K, which is not apparent. Alternatively, the cold dust may be located 
in a very hot but diffuse phase of the gas, surrounding the bright, dense X-ray filaments. If the gas 
had a density of < 0.01 cm~^ then the grains could be this cold. The dust could also be contained 
in dense clumps which are in not in pressure equilibrium with the more diffuse X-ray gas, a model 
which has been suggested for SN 1987A^^. Finally, the cold dust could be in thermal equilibrium 
with the same X-ray gas which is responsible for heating the hot dust, if the cold grains are both 
very large (1-10/xm) and emit more efficiently than grains in the diffuse ISM. Evidence for unusual 
grain properties in Cas A is presented below. 

The mass of dust can be estimated from the submm emission using = ^ ^b{vT) ' where D is the 
distance to the remnant (3.4 kpc), B{u, T) is the Planck function at temperature T and k is the dust 
mass absorption coefficient. The submm emission is clearly dominated by the cold dust, therefore the 
cold temperature is the only one we need be concerned with. The value of k is the main uncertainty 
as its value is not well determined and may vary with physical environment. Dust masses derived 
from the SED and different values of n are in Table 1. Using 'standard' k, values appropriate to the 
diffuse ISM^*^ [n^) gives dust masses which are uncomfortably high (> TM©). Higher n values, as 
measured in environments where dust may be newly formed, amorphous or coagulated (see Table 
1) produce more sensible estimates (2-6 M©), thus the dust in Cas A appears to have different 
properties to that in the diffuse ISM. It is possible that grain processing in the diffuse ISM by UV 
photons and passage through dark clouds may alter the surface chemistry or shape of grains, causing 
a change in their cmissivity. At an age of 320 yrs, the dust in Cas A is relatively pristine and has 
not yet been subjected to the above events, this may contribute to its greater emissivity. 

Was the dust produced in the supernova explosion, or could it be pre-existing dust which was 
swept up by the blast wave? If the remnant has swept up ISM material with a uniform density^^ of 
0.4-4 H atoms per cm~^ then only 0.004 — 0.04 Mq of dust should be present^ Higher ISM densities 
are possible, but in order for the mass of dust inferred in Cas A to have been swept up densities of the 
order 200-400 cm^^ are required. This would mean Cas A has swept up 200 — 400 Mq of material 
which is inconsistent with its dynamical state^^'^^, believed to be just entering the Sedov stage (i.e. 
it has swept up about as much mass as it ejected (< ISM©)). Another hypothesis for the bright ring 
seen in the X-ray and radio is that it is the interaction of the supernova with a stellar wind shell 
swept up by the Wolf-Rayet precursor In this case, one could estimate that somewhere between 
4-16 Mq of material could have been lost in the red-supergiant phase and swept up into the ring. 
For this to be the source of the dust we observe requires an implausibly large quantity of freshly 

synthesised elements efficiently condensing into dust in the stellar wind material^ (which, unlike the 
Hot an ISM gas-to-dust ratio of 160 
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supernova ejecta, is mainly H and He). We conclude that the most likely explanation for the large 
mass of dust in Cas A is that it was formed as a result of the supernova explosion. 

Using our result, we estimate the current Galactic dust production rate from Typc-II SNc as 
7 — 18 (x 10~^) M0 yr^-*^. Stellar wind^ sources produce ~ 5 x lO~^M0yr~^, meaning that SNc are 
the dominant source of interstellar dust during most of galactic evolution. We have shown here that 
individual SNe are capable of producing of order of a solar mass of dust in the short times available 
in the early universe and therefore the dust observed at high redshift is likely to have originated in 
super novae. 
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Table 1: Properties of Cas A 

SED parameters Dust masses (Mq) 

P Thot (K) Tcold (K) Ki K2 K3 

0.91°^ 11211} IS.Otli M45o = 2.2li| M450 = 3.81^° 71^450 = 12.71?:? 

M85o = 2.4l°:^ Ms5o = 6.2tli Ms5o = 26.0tli 

SED parameters are the best least-squares fit to the data points. The errors were estimated using a 
bootstrap procedure in which the measured fluxes were perturbed according to their Gaussian 
errors and re-fitted to find new parameters. The process was repeated 3000 times and the 
uncertainties quoted are the 68% confidence interval derived from this method. Quoted dust 
masses are those from the best-fit SED parameters for each value of k, and the ranges refiect the 
68% confidence interval using the distribution produced by the bootstrap, essentially the 
uncertainty in dust mass due to the uncertainty in fluxes and SED fltting. We have used three 
values of k which represent the average from numerous source in the literature for different grain 
environments (to be presented in more detail elsewhere, Dunne et al. in prep). (1) Laboratory 
measurements and theoretical modelling suggest that if grains are amorphous or have a clumpy, 
aggregate structure k can be very high 0.6 — l.lm^kg^^ at 850fim. We take 
Ki(450) = 1.5m^kg~i, ki(850) = 0.76 m^kg""^. (2) Environments where dust may be newly 
formed, or modified by icy mantles also show high values 0.16 — 0.8m^kg~i. From observations of 
evolved stars, VRN (visual reflection nebulae), planetary nebulae and cold clouds in the Galaxy we 
take K2(450) = 0.88 m^ kg"^ , «;2(850) = 0.3 m^ kg"^. (3) Studies of extragalactic systems and 
diffuse ISM dust in the Galaxy support lower values of k in the region of 0.04 — 0.15 m^kg"^ at 
850/xm. We take ^3(450) = 0.26 m^ kg-^ , ^3(860) = 0.073 m^ kg-^. 
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Figure 1: SCUBA 850//m image of Cas A at a resolution of 15 arcsec. The data were taken in scan 
map mode and reduced in a standard way using SURF. Scan map is known to produce artifacts on 
large angular scales on the final image and these remain the limitation on the accuracy of the absolute 
flux measurement and morphology. On the 850/im image, the linear baseline removal left an area 
of diffuse emission to the west of the remnant which we believe was an artifact. After subtraction 
of a surface which removed this feature the loss of flux in the remnant was 5 Jy. We confirmed 
the accuracy of the absolute flux level of the maps by making new photometric measurements with 
SCUBA at several positions on the remnant in December 2002. The quoted la uncertainty on the 
fluxes is a combination of calibration errors (5-10% and 15-20% at 850/450/im), plus the range of 
fluxes from the maps which agreed with the photometry. The integrated fluxes in were corrected for 
the error lobe contribution (derived from maps of Uranus) by dividing by 1.2 at 850/xm and 2.0 at 
450jUm. 
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Figure 2: The SED of Cas A from the mid-IR to the radio. All fluxes are integrated. The radio 
spectrum is fitted with a power law slope with spectral index a = —0.72. The radio data^*^'^"^^^^ have 
been scaled to the epoch of our Cas A observation (1998.4) using a recent study^^ which found that 
the flux decrease of Cas A was independent of frequency at 0.6 — 0.7%yr~^. The IRAS points are 
the averages of the literature values^^'^^ with the error bars indicating the range of fluxes. SCUBA 
fluxes are 50.8 ± 5.6 Jy at 850/xm and 69.8 it 16.1 Jy at 450/xm. The synchrotron contribution is 
34.9 Jy at 850//m and 22.1 Jy at 450//m. The fit to the FIR/submm points is for a 2 temperature 
grey-body with (3 = 0.9, T^, = 112 K, = 18 K and for 700 times more mass in the cold component 
than in the hot. We have fitted 5 parameters to 6 data points and therefore the fit is just constrained. 
The uncertainties in the SED parameters are given in Table 1. Clearly, the submm points at 450 
and 850/Ltm lie above the extrapolation of the radio synchrotron spectrum. The point at 1.2mm (ref. 
28) is not of sufficient accuracy to determine an excess, but is consistent with our higher frequency 
measurements. Longer wavelength ISO measurements^'^ of Cas A at 170/xm suggested that dust 
at ^ 30 K may have been present (implying of order 0.15 Mq of dust), however the poor angular 
resolution of ISO at these wavelengths made a separation of the remnant and background emission 
impossible. The ISO fiux is therefore a lower limit (R. J. Tuffs, private communication). 
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Figure 3: SCUBA 450/im map, smoothed with a 21" gaussian. Colours and contours both represent 
the 450//m emission. Contours start at 2a, with intervals of 1 a. The difference in morphology 
between this image and that at 850/um (Fig. 1) is because emission from cold dust dominates at 
450//m, while two-thirds of the emission at 850/im is synchrotron. 
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Figure 4: The 850/um emission once the synchrotron has been subtracted using an 83 GHz image . 
The box is 8.4 arcmin by 7.8 arcmin, north is up and east is left. Colours represent the 850/im 
intensity, contours are the 450;um emission with the synchrotron subtracted, starting at 3cr with 
increments of +lcr. The rings and crosses indicate the location (and centroids) of the forward and 
reverse shocks as determined from Chandra X-ray data^^. The forward shock is at a mean radius of 
153 lb 12" and the reverse shock is at a mean radius of 95± 10". The bulk of the dust emission appears 
to be bounded by the shocks, where the gas density is highest. There is a noticeable asymmetry in 
the dust emission. 
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